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The binding of a ligand molecule to a protein is often accompanied by conformational changes of the protein. 
A central question is whether the ligand induces the conformational change (induced-fit), or rather selects and 
stabilizes a complementary conformation from a pre-existing equilibrium of ground and excited states of the 
protein (selected-fit). We consider here the binding kinetics in a simple four-state model of ligand-protein bind- 
ing. In this model, the protein has two conformations, which can both bind the ligand. The first conformation is 
the ground state of the protein when the ligand is off, and the second conformation is the ground state when the 
ligand is bound. The induced-fit mechanism corresponds to ligand binding in the unbound ground state, and the 
selected-fit mechanism to ligand binding in the excited state. We find a simple, characteristic difference between 
the on- and off-rates in the two mechanisms if the conformational relaxation into the ground states is fast. In the 
case of selected-fit binding, the on-rate depends on the conformational equilibrium constant, while the off-rate 
is independent. In the case of induced-fit binding, in contrast, the off-rate depends on the conformational equi- 
librium, while the on-rate is independent. Whether a protein binds a ligand via selected-fit or induced-fit thus 
may be revealed by mutations far from the protein's binding pocket, or other "perturbations" that only affect the 
conformational equilibrium. In the case of selected-fit, such mutations will only change the on-rate, and in the 
case of induced-fit, only the off-rate. 



INTRODUCTION 

The function of a protein often requires dynamic changes 
between different three-dimensional conformations. Confor- 
mational changes that occur when a protein binds a ligand 
molecule are well documented by experimentally determined 
structures of proteins in the unbound state and the ligand- 
bound state illEl. In 1958, Koshland |3| suggested that the 
binding of the ligand may induce the conformational change 
of the protein. More recently, Tsai et al. [4,, JJ suggested an 
alternative mechanism in which the ligand selects and stabi- 
lizes a complementary protein conformation from an equilib- 
rium of low-energy and higher-energy conformations. This 
selected-fit, or conformational selection mechanism is based 
on the energy-landscape picture of protein folding 

How can we distinguish whether a protein binds its ligand 
in an induced-fit or selected-fit mechanism? Recent single- 
molecule fluorescence lH HI [TOl [IT] [ill [IS and NMR relax- 
ation experiments lfl4l [TSl [T6l [TTl [TSl provide dynamic in- 
sights into conformational transitions. NMR relaxation exper- 
iments of the enzyme DHFR, which binds two ligands, reveal 
excited-state conformations that resemble ground-state con- 
formations of adjacent ligand-bound or ligand-free states on 
the catalytic cycle |16|. These experiments provide evidence 
for a pre-existing equilibrium of conformations that is shifted 
by the binding or unbinding of a ligand, as suggested in the 
selected-fit mechanism of protein binding H^l. For the bind- 
ing of myosin to actin, in contrast, an induced-fit mechanism 
has been proposed fT\ based on cryo-microscopic images of 
the actin-myosin complex |[T9l . 

We consider here the selected-fit and induced-fit binding 
kinetics in a simple four-state model of protein-ligand bind- 
ing (see fig. [T}. In this model, the protein, or enzyme, has 
two dominant conformations Ei and £2- The conformation 
El is the ground-state conformation in the unbound state of 
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FIG. 1 : Simple four-state model of protein-ligand binding. Without 
the ligand L, the conformation E\ of the protein, or enzyme, is the 
ground-state conformation, and E2 is the excited-state conformation. 
When the ligand is bound, E2L is the ground state, and E\L the 
excited state. On the selected-fit route, the ligand binds the protein 
in the excited-state conformation E2, and on the induced-fit binding 
route, in the ground-state conformation E\ . 



the protein, while E2 is the ground-state conformation in the 
ligand-bound state. Two routes connect the unbound ground 
state El and bound ground state E2L. On the induced-fit route 
El EiL ^ E2L, the protein first binds the ligand in con- 
formation El, which causes the transition into conformation 
E2. On the selected-fit route Ei E2 E2L, the protein 
binds the ligand in the higher-energy conformation E2- 

We find a characteristic difference between the selected-fit 
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and induced-fit binding kinetics. If the conformational relax- 
ation into the ground state is fast, the selected-fit on-rate de- 
pends on the equilibrium constant of the conformations Ei 
and E2, while the selected-fit off-rate is independent of the 
conformational equilibrium. The induced-fit on-rate, in con- 
trast, is independent of the conformational equilibrium be- 
tween EiL and E2L, whereas the induced-fit off -rate depends 
on this equilibrium. Mutations or other perturbations that shift 
the conformational equilibrium without affecting the shape 
or free energies of the binding site thus may help to identify 
whether a protein binds its ligand in a selected-fit or induced- 
fit mechanism. 



RESULTS 

Equilibrium behavior 

The four states of our model are connected by four tran- 
sitions with equilibrium constants Ki, K2, K^, and Ki, (see 
fig. [T}- In the model, the constants and Kb for the con- 
formational equilibrium in the unbound and the bound state 
obey 



[El] 



< 1 



[E2L] 
[EiL] 



> 1 



(1) 



since Ei is the ground-state conformation in the unbound state 
of the protein, and E2L is the ground state when the ligand is 
bound. We assume here that the excited-state conformations 
are significantly higher in free energy than the ground states. 
The binding equilibrium of the two conformations is governed 
by the constants 



[EiL] 
[Ei][L] 



[E2L] 

[E2m 



(2) 



From the two equahties [E2L] — K2[E2][L] = 
K2Ku[Ei][L] and [E2L] = Kb[EiL] = KbKi[Ei][L] that 
follow directly from these definitions, we obtain the general 
relation 



KuK2 = KbKi 



(3) 



between the four equilibrium constants. Thus, only three of 
the equilibrium constants are independent. 

The selected-fit route Ei E2 E2L in our model 
dominates over the induced-fit route Ei EiL E2L if 
the concentration [E2] of the selected-fit intermediate is larger 
than the concentration [EiL] of the induced-fit intermediate. 
Since [E2] = K^[Ei] and [EiL] = Ki[Ei][L], the selected- 
fit mechanism thus is dominant for small ligand concentra- 
tions 



[L] < K^/K, 



(4) 



while the induced-fit mechanism is dominant for large ligand 
concentrations [L] > Ku/Ki. 



Selected-fit binding kinetics 

We first consider the binding kinetics along the selected-fit 
route of our model (see figs.[T]and[2]). The selected-fit binding 
rate is the dominant relaxation rate of the process 



El 



st[L] 
E2 ^ E2L 



(5) 



from El to E2L, with E2L as an "absorbing state" with- 
out backflow into i?2- Here, S21 is the transition rate from 
state El to E2, S12 is the rate for the reverse transition, and 
Sb is the binding rate per mole ligand in state E2- In the 
appendix, we calculate the exact relaxation rates for a pro- 
cess of the form (|5]l. Since E2 is the excited state, we have 
S21/S12 = Ku <C 1, and the on-rate of the selected-fit pro- 
cess (|5]) is approximately given by 



S2iSb[E] 
ML] 



S12 



(6) 



(see eq. p5| ) in the Appendix). For small ligand concen- 
trations [L], or fast conformational relaxation S[2 into the 
ground-state conformation Ei, we have Sb[L] ^ sn, and the 
selected-fit on-rate per mole ligand is 



Sl2 



(7) 



The selected-fit on-rate (j?]) just depends on the equilibrium 
constant of the conformations Ei and E2 in the unbound 
state, and the binding rate Sb of the conformation E2- Since 
the equilibrium probabiUty P{E2) = [E2]/{[Ei] + [E2]) of 
conformation is approximately P(£^2) ~ [^'2]/[-£'i] = Ku 
for [E2] ^ [-El], the selected-fit on-rate ^ can also be di- 
rectly understood as the product of the probabiUty that the 
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FIG. 2: Binding kinetics along the selected-fit route of our model 
(see fig. [T}. Here, S21 and S12 are the rates for the conformational 
transitions in the unbound state, sj, is the binding rate of conforma- 
tion E2 per mole ligand, and s,, the unbinding rate. Since Ei is the 
ground state and E2 the excited state, we have Ku ~ S21/S12 <^ 1. 
If the conformational transition rate S12 into the ground state is much 
larger than the binding and unbinding rates st[L] and s„, the on-rate 
along the selected-fit route is approximately Son ~ K^Sb, and the 
off-rate is Soff ~ s„ (see eqs. |7]l and |9]l). 
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protein is in conformation E2 and the binding rate Sf, of this 
conformation. 

The selected-fit off -rate is the dominant relaxation rate of 
the process 



' on 



Su Sl2 

EoL -r — ^ Eo ^ E-i 



(8) 



from E2L to the Ei. This process follows the same general 
reaction scheme as the binding process (j5|. A reasonable, sim- 
plifying assumption here is that the conformational relaxation 
process from the excited state E2 to the ground state Ei is 
significantly faster than the binding and unbinding process, 
which implies S12 3> and S12 ^ S}j[L\. The selected-fit 
off-rate then simply is 



Soff ' 



(9) 



(see eq. ( 26 1 in the Appendix). The off -rate Soff thus is inde- 
pendent of the conformational transition rates between Ei and 
i?2. The off -rate is identical with the rate for the bottleneck 
step, the unbinding process from E2L to i?2. 



Induced-fit binding kinetics 

The on-rate along the induced-fit binding route of our 
model (see figs. [T] and [3]l is the dominant relaxation rate of 
the process 



rb[L ] r2i ^ 

El — ^ E\L ^ E2L 



(10) 



Here, r^, is the binding rate of conformation Ei per mole lig- 
and, r„ is the unbinding rate, and r2i is the rate for the con- 
formational transition into the bound ground state E2L. The 



induced-fit binding process ( lOi is similar to the selected-fit 
unbinding process (|8]l. As before, we assume that the con- 
formational transition into the ground state is much faster 
than the binding and unbinding processes, i.e. we assume 
^21^ rb[L\ and r2i ^ r^- The induced-fit on-rate per mole 
ligand then is (see Appendix) 



(11) 



The induced-fit unbinding rate is the dominant relaxation 
rate of the process 



1'12 Tu 

E2L ^ — ^ EiL ^ El 

r2i 



(12) 



which is similar to the selected-fit binding process (|5]). Since 
E2L is the ground-state conformation in the bound state of 
the protein, we have Kb — r2\jrYi ^ 1- The dominant re- 
laxation rate of ([12]) then is 



^off ' 



r2i + rt 



(13) 



> rb[l-] ^ 



'21 



12 




'off 



FIG. 3: Binding kinetics along the induced-fit route of our model 
(see fig.[T](. Here, ri2 and r2i are the rates for the conformational 
transitions in the bound state, rt is the binding rate of conformation 
El per mole ligand, and the unbinding rate. Because E2L is the 
ground state, we have Kb = r2i/rr2 ^ 1- For conformational 
transition rates 7-21 into the bound ground state that are much larger 
than the binding and unbinding rates rb[L] and r„, the on- and off- 
rates along the induced-fit route are approximately ron ~ rt, and 
'"off ~ Sii/ Kb (see eqs. 1 1 1 1 and 1 14 1). 



(see eq. (25i in the Appendix). For fast conformational re- 
laxation into the bound ground state with rate r2i ^ r^, the 
induced-fit off -rate is approximately 



''off ■ 



Kb 



(14) 



The off -rate r^jKb can again be understood as the product of 
the unbinding rate in the excited state E^L and the equilib- 
rium probability P{EiL) = [EiL]/{[EiL] + [E2L]) of this 
state, which is approximately P{EiL) « [EiL]/[E2L] = 
l/Kb. 



Mutational analysis of the binding kinetics 

A mutational analysis of the binding kinetics may reveal 
the characteristic differences between the selected-fit and 
induced-fit on-rates (|7| and ( [TT| and off-rates (|9]) and ( [T4] l. 
Of particular interest here are mutations of residues far away 
from the ligand-binding site that affect the rate constants for 
transitions between the conformations, but not the binding 
rates Sb and Vb and unbinding rates Su and r„ of the two con- 
formations. The mutations will change the free-energy dif- 
ferences AG„ = G{E2) - G{Ei) and AGb = G{E2L) - 
G{EiL) between the two conformations in the unbound and 
the bound state. In terms of these free-energy differences, the 
equilibrium constants Ku and Kb for the conformational tran- 
sitions are 

= i?rexp(-AG„/i?T) and K^ = RTexp{^AGb/RT) 

(15) 

From eq. (|7]i, we find that the ratio of the selected-fit on-rates 
for wildtype and mutant 



— = §^= cxp{-AAG JRT) 



(16) 
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depends on the mutation-induced shift AAG^ = AG^ — 
AGtj of the free energy difference between the conformations 
in the unbound state. The prime here indicates the conforma- 
tional free-energy difference for the mutant. From eq. ( 14 1, 
we obtain the relation 



'off 
foff 



Kb 
KL 



exp(AAGb/i?T) 



(17) 



between the induced-fit off-rates of wildtype and mutant. 
Here, AAG;, — AG[, — AG^ is the mutation-induced shift of 
the conformational free-energy difference in the bound state. 
The selected-fit off-rate (|9]) and induced-fit on-rate ( [TT| , in 
contrast, are the same for wildtype and mutant if the muta- 
tions do not affect the binding and unbinding rates of the two 
conformations. 



DISCUSSION AND CONCLUSIONS 

Empirical methods to calculate mutation-induced stability 
changes of proteins have been investigated intensively in the 
pastyeai-s lIS|2ll|22l|23l|24l|25l|2g|^|28l|29| These meth- 
ods can be used to calculate the changes AAG„ of the con- 
formational free-energy difference in the unbound state, pro- 
vided that experimental protein structures are available both 
for the unbound and ligand-bound state of the protein. The 
experimental structure of the protein in the ligand-bound state 
then corresponds to the excited-state conformation in the un- 
bound state of our model. Calculations of mutation-induced 
changes AAGb of conformational free-energy differences in 
the bound state are more difficult and require a construction 
of a ligand-bound excited-state conformation by docking the 
ligand to the experimental structure of the unbound protein. 

In combination with eqs. ( 16 1 and (17 1 of our model, cal- 



culated free-energy changes may help to identify selected-fit 
or induced-fit mechanisms of protein binding. A promising 
candidate is the protein DHFR, which has been studied exten- 
sively as a model enzyme to understand the relations between 
conformational dynamics, binding, and function. The enzy- 
matic mechanism and conformational changes of DHFR have 
been investigated by mutational analyses ll30ll3Tll32ll33l[34l . 
NMR relaxation |16, 17 35 1, single-molecule fluorescence 
experiments |10|, and simulation i36ll37ll38l . The strong ef- 
fect of point mutations far from the ligand binding site on cat- 
alytic rates and binding rates suggests large conformational 
changes during the catalytic cycle |[30l [32l l33l [34l |39l , in 
agreement with experimental structures of DHFR in the un- 
bound and several ligand-bound states ll40l . Mutations far 
from the binding sites of DHFR affect the binding kinetics of 
the two ligands 113111321 |33J . and NMR relaxation experiments 
point towards a selected-fit mechanism of ligand binding lfT6l . 
However, the mutations seem to affect both the binding and 
the unbinding rates |31 , 32. 33 1. This precludes a direct com- 
parison with our model, either because the conformational dy- 
namics is more complex, or because the mutations have an 
indirect effect on the binding free energies. 

We have considered here a simple four-state model of 
protein-ligand binding. The kinetic and mutational analysis 



of this model is partly inspired by previous models of the 
protein folding kinetics BTl |42 I. We have found character- 
istic kinetic differences between the selected-fit and induced- 
fit binding routes of this model, which should be observ- 
able also in more detailed models of conformational changes 
and binding. While the timescale of the relevant confor- 
mational transitions are in general inaccessible to standard 
molecular dynamics simulations with atomistic protein mod- 
els |43|, conformational selection during binding has been re- 
cently studied by combining atomistic simulations and dock- 
ing fl4l [45] |46l . Large conformational transitions are also 
intensively studied by normal mode analysis |47, 48 1, in par- 
ticular with coarse-grained elastic network models of proteins 

lllSlEolEIlEIlEaEai. 



APPENDIX: RELAXATION RATES 

General solution 

The binding and unbinding processes of our model have a 
general reaction scheme of the form 



B 



C 



(18) 



For the selected-fit or induced-fit binding process, state A cor- 
responds to El, and state G to E2L. For the unbinding pro- 
cesses, A corresponds to E2L, and G to £'1. We are interested 
in the dominant, slowest relaxation rate into state G. There- 
fore, the state G is an "absorbing" state, i.e. the rate from G 
back to B is 0. The probability evolution PA{t), ^s(i), and 
Pc{t) of the three states is then governed by the set of master 
equations 



dPAjt) 
dt 

dt 
dPcjt) 
dt 



= kABPB{t) ~ kBAPA{t) (19) 
= kBAPA{t)-{kAB + kcB)PB{t) (20) 
= kcBPBit) (21) 



where kij denotes the rate from state j to state i. The three 
equations can be written in the matrix form 



dPjt) 
dt 



-WP{t) 



with P{t) = {PA{t),PB{t),Pc{t)) and 

kBA -kAB 0^ 

W = I -kBA kAB + kcB 

-kcB Oj 



The general solution of ( 22 1 has the form 115511 

Pit) = CoYo + ciYie-^i* + C2Yie-^^* 



(22) 



(23) 



(24) 



where Yo, Yi, and Y2 are the three eigenvectors of the ma- 
trix W, and Ai and A2 are the two positive eigenvalues of 
the eigenvectors Yi and Y2- The coefficients Co, ci, and C2 
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depend on the initial conditions at time t = 0. The eigen- 
value of the eigenvector Yo is 0, which ensures that P{t) 
relaxes towards a finite equilibrium probability distribution 
P{t oo) — CqYo for t ^ oo. The two nonzero eigen- 
values are 



1 



± 



{ksA - kcBY + 2kAB{kBA + kcB) 



These eigenvalues, and their eigenvectors Yi and Y2, capture 
the time-dependent relaxation process into equilibrium (see 
eq. my 



P(0) = (1,0,0) (only state A populated) or P(0) = (1 - 
kBA/kAB,kBA/kAB,0) ('pre-equilibration' of states A and 
B). For both initial conditions, the probability evolution of 
state C is approximately given by Pc{t) ~ 1 — cxp[— Ait]. 
The dominant relaxation rate Ai of eq. ( p5] l can also be ob- 
tained from eqs. ( [T9| to pT| ) in a steady-state approximation, 
i.e. under the assumption that the variation dPB{t)/dt of the 
intermediate state B is negligibly small. 



Selected-fit unbinding and induced-fit binding 



Selected-fit binding and induced-fit unbinding 

In selected-fit binding (|5]l, state A corresponds to Ei, state 
B to E2, and C to E2L (see fig.[T]l. In induced-fit unbinding 
(fT2]i, A coiTesponds to E2L, B to EiL and C to Ei. In our 
model, we have [E2] ^ [Ei] and [E2L] :$> [EiL] in equilib- 
rium, which implies kg^ kAB in both cases. The nonzero 
eigenvalues then are approximately 



Ai 



kBAkcB 
kAB + kcB 



and A 9 



kAB 



^CB 



(25) 



For kBA ^ kAB, we have Ai <C A2. The smaller rate 
Ai is the dominant relaxation rate for the initial conditions 



In selected-fit unbinding (|8]l, the rate kcB corresponds to 
S12, the rate from the excited state E2 in the unbound ground 



state El. In induced-fit binding ( 10 1, kcB corresponds to r2i, 
the rate from EiL in the bound ground state E2L. We assume 
here that the rates for these conformational transitions into the 
ground state are much larger than the binding and unbinding 
rates. This implies kAB <C fees and ksA <C fees- The 
nonzero eigenvalues then simply are 



Ai w kBA and A2 « fees 



(26) 



with the clearly smaller eigenvalue Ai as the dominant rate of 
the relaxation process from state A to C. 



[1] Gerstein, M. and Krebs, W. A database of macromolecular mo- 
tions. Nucleic Acids Res. 26:4280-4290, 1998. 

[2] Goh, C.-S., Milbum, D., and Gerstein, M. Conformational 
changes associated with protein-protein interactions. Curr. 
Opin. Struct. Biol. 14:104-109, 2004. 

[3] Koshland, D. E. Application of a theory of enzyme specificity 
to protein synthesis. Proc. Natl. Acad. Sci. USA 44:98-104, 
1958. 

[4] Tsai, C. J., Kumar, S., Ma, B., and Nussinov, R. Folding 
funnels, binding funnels, and protein function. Protein Sci. 
8:1181-1190, 1999. 

[5] Tsai, C. J., Ma, B., Sham, Y. Y., Kumar, S., and Nussinov, 
R. Structured disorder and conformational selection. Proteins 
44:418-427,2001. 

[6] Dill, K. A. and Chan, H. S. From Levinthal to pathways to 
funnels. Nat. Struct. Biol. 4:10-19, 1997. 

[7] Bryngelson, J. D., Onuchic, J. N., Socci, N. D., and Wolynes, 
P. G. Funnels, pathways, and the energy landscape of protein 
folding: a synthesis. Proteins 21:167-195, 1995. 

[8] Lu, H. P., Xun, L., and Xie, X. S. Single-molecule enzymatic 
dynamics. Science 282: 1877-1882, 1998. 

[9] Ha, T, Ting, A. Y, Liang, J., Caldwell, W. B., Deniz, A. A., 
Chemla, D. S., Schultz, P. G., and Weiss, S. Single-molecule 
fluorescence spectroscopy of enzyme conformational dynamics 
and cleavage mechanism. Proc. Natl. Acad. Sci. USA 96:893- 
898, 1999. 

[10] Antikainen, N. M., Smiley, R. D., Benkovic, S. J., and Hammes, 



G. G. Conformation coupled enzyme catalysis: single-molecule 

and transient kinetics investigation of dihydrofolate reductase. 

Biochemistry 44:16835-16843, 2005. 
[11] Min, W., EngUsh, B. P, Luo, G., Cherayil, B. J., Kou, S. C, and 

Xie, X. S. Fluctuating enzymes: lessons from single-molecule 

studies. Acc. Chem. Res. 38:923-931, 2005. 
[12] Michalet, X., Weiss, S., and Jager, M. Single-molecule fluores- 
cence studies of protein folding and conformational dynamics. 

Chem. Rev. 106:1785-1813, 2006. 
[13] Smiley, R. D. and Hammes, G. G. Single molecule studies of 

enzyme mechanisms. Chem. Rev. 106:3080-3094,2006. 
[14] Eisenmesser, E. Z., Bosco, D. A., Akke, M., and Kern, D. 

Enzyme dynamics during catalysis. Science 295:1520-1523, 

2002. 

[15] Wang, C, Karpowich, N., Hunt, J. F, Ranee, M., and Palmer, 
A. G. Dynamics of ATP-binding cassette contribute to allosteric 
control, nucleotide binding and energy transduction in ABC 
transporters. J. Mol. Biol. 342:525-537, 2004. 

[16] Boehr, D. D., McElheny, D., Dyson, H. J., and Wright, P E. The 
dynamic energy landscape of dihydrofolate reductase catalysis. 
Science 313:1638-1642, 2006. 

[17] Boehr, D. D., Dyson, H. J., and Wright, P E. An NMR perspec- 
tive on enzyme dynamics. Chem. Rev. 106:3055-3079, 2006. 

[18] Henzler-Wildman, K. and Kern, D. Dynamic personalities of 
proteins. Nature 450:964-972, 2007. 

[19] Holmes, K. C, Angert, I., KuU, F J., Jahn, W., and Schroder, 
R. R. Electron cryo-microscopy shows how strong binding of 



6 



myosin to actin releases nucleotide. Nature 425:423^27, 2003. 

[20] Guerois, R., Nielsen, J. E., and Serrano, L. Predicting changes 
in the stability of proteins and protein complexes: a study of 
more than 1000 mutations. J. Mol. Biol. 320:369-387, 2002. 

[21] Schymkowitz, J., Borg, J., Stiicher, R, Nys, R., Rousseau, R, 
and Serrano, L. The FoldX web server: an online force field. 
Nucleic Acids Res. 33:W382-W388, 2005. 

[22] Gromiha, M. M., Oobatake, M., Kono, H., Uedaira, H., and 
Sarai, A. Role of structural and sequence information in the pre- 
diction of protein stability changes: comparison between buried 
and partially buried mutations. Protein Eng. 12:549-555, 1999. 

[23] Gilis, D. and Rooman, M. PoPMuSiC, an algorithm for pre- 
dicting protein mutant stability changes: application to prion 
proteins. Protein Eng. 13:849-856, 2000. 

[24] Zhou, H. and Zhou, Y. Distance-scaled, finite ideal-gas refer- 
ence state improves structure-derived potentials of mean force 
for structure selection and stability prediction. Protein. Sci. 
11:2714-2726, 2002. 

[25] Capriotti, E., Rariselli, R, and Casadio, R. I-Mutant2.0: predict- 
ing stability changes upon mutation from the protein sequence 
or structure. Nucleic Acids Res. 33:W306-W310, 2005. 

[26] Cheng, J., Randall, A., and Baldi, P. Prediction of protein stabil- 
ity changes for single-site mutations using support vector ma- 
chines. Proteins 62:1125-1132, 2006. 

[27] Huang, L.-T., Saraboji, K., Ho, S.-Y., Hwang, S.-R, Pon- 
nuswamy, M. N., and Gromiha, M. M. Prediction of protein 
mutant stability using classification and regression tool. Bio- 
phys. Chem. 125:462^70, 2007. 

[28] Parthiban, V., Gromiha, M. M., Hoppe, C., and Schomburg, D. 
Structural analysis and prediction of protein mutant stability us- 
ing distance and torsion potentials: role of secondary structure 
and solvent accessibility. Proteins 66:41-52, 2007. 

[29] Bueno, M., Camacho, C, and Sancho, J. SIMPLE estimate of 
the free energy change due to aliphatic mutations: Superior pre- 
dictions based on first principles. Proteins 68:850-862, 2007. 

[30] Cameron, C. E. and Benkovic, S. J. Evidence for a functional 
role of the dynamics of glycine- 121 of Escherichia coli dihydro- 
folate reductase obtained from kinetic analysis of a site-directed 
mutant. Biochemistry 36: 15792-15800, 1997. 

[31] Miller, G. P. and Benkovic, S. J. Deletion of a highly motional 
residue affects formation of the Michaelis complex for Es- 
cherichia coli dihydrofolate reductase. Biocherrustry 37:6327- 
6335, 1998. 

[32] Miller, G. P. and Benkovic, S. J. Strength of an interloop hy- 
drogen bond determines the kinetic pathway in catalysis by Es- 
cherichia coU dihydrofolate reductase. Biochemistry 37:6336- 
6342, 1998. 

[33] Miller, G. P., Wahnon, D. C, and Benkovic, S. J. Interloop con- 
tacts modulate ligand cycling during catalysis by Escherichia 
coli dihydrofolate reductase. Biochemistry 40:867-875, 2001. 

[34] Wang, L., Goodey, N. M., Benkovic, S. J., and Kohen, A. Co- 
ordinated effects of distal mutations on environmentally cou- 
pled tunneling in dihydrofolate reductase. Proc. Natl. Acad. 
Sci. USA 103:15753-15758, 2006. 

[35] McElheny, D., Schnell, J. R., Lansing, J. C, Dyson, H. J., and 
Wright, P. E. Defining the role of active-site loop fluctuations in 
dihydrofolate reductase catalysis. Proc. Natl. Acad. Sci. USA 
102(14):5032-5037, 2005. 

[36] Agarwal, R K., Billeter, S. R., Rajagopalan, R T. R., Benkovic, 
S. J., and Hammes-Schiffer, S. Network of coupled promot- 
ing motions in enzyme catalysis. Proc. Natl. Acad. Sci. USA 
99:2794-2799, 2002. 

[37] Rod, T. H., Radkiewicz, J. L., and Brooks III, C. L. Corre- 
lated motion and the effect of distal mutations in dihydrofolate 



reductase. Proc. Nafl. Acad. Sci. USA 100:6980-6985, 2003. 

[38] Chen, J., Dima, R. I., and Thirumalai, D. AUosteric communi- 
cation in dihydrofolate reductase: signaling network and path- 
ways for closed to occluded transition and back. J. Mol. Biol. 
374:250-266, 2007. 

[39] Benkovic, S. J., Hammes, G. G., and Hammes-Schiffer, S. Free- 
energy landscape of enzyme catalysis. Biochemistry 47:3317- 
3321, 2008. 

[40] Sawaya, M. R. and Kraut, J. Loop and subdomain movements 
in the mechanism of Escherichia coli dihydrofolate reductase: 
crystallographic evidence. Biochemistry 36:586-603, 1997. 

[41] Merlo, C, Dill, K. A., and Weikl, T. R. ^-values in protein fold- 
ing kinetics have energetic and structural components. Proc. 
Nafl. Acad. Sci. USA 102:10171-10175, 2005. 

[42] Weikl, T. R. Transition states in protein folding kinetics: Mod- 
eling $-values of small /3-sheet proteins. Biophys. J. 94:929- 
937, 2008. 

[43] Karplus, M. and McCairmion, J. A. Molecular dynamics simu- 
lations of biomolecules. Nat. Struct. Biol. 9:646-652, 2002. 

[44] Wong, S. and Jacobson, M. P. Conformational selection in sil- 
ico: loop latching motions and ligand binding in enzymes. Pro- 
teins 71:153-164, 2008. 

[45] Wong, C. R, Kua, J., Zhang, Y, Straatsma, T. R, and McCam- 
mon, J. A. Molecular docking of balanol to dynamics snapshots 
of protein kinase A. Proteins 61:850-858, 2005. 

[46] Frembgen-Kesner, T. and Elcock, A. H. Computational sam- 
pling of a cryptic drug binding site in a protein receptor: ex- 
plicit solvent molecular dynamics and inhibitor docking to p38 
MAP kinase. J. Mol. Biol. 359:202-214, 2006. 

[47] Ma, J. Usefulness and limitations of normal mode analysis 
in modeling dynamics of biomolecular complexes. Structure 
13:373-380, 2005. 

[48] Case, D. Normal mode analysis of protein dynamics. Curr. 
Opin. Struct. Biol. 4:285-290, 1994. 

[49] Bahar, I. and Rader, A. J. Coarse-grained normal mode analy- 
sis in structural biology. Curr. Opin. Struct. Biol. 15:586-592, 
2005. 

[50] Yang, L., Song, G., and Jemigan, R. L. How well can we imder- 
stand large-scale protein motions using normal modes of elastic 
network models? Biophys. J. 93:920-929, 2007. 

[51] Tirion, M. Large amplitude elastic motions in proteins from a 
single-parameter, atomic analysis. Rhys. Rev. Lett. 77:1905- 
1908, 1996. 

[52] Bahar, I., Atilgan, A. R., and Erman, B. Direct evaluation of 
thermal fluctuations in proteins using a single-parameter har- 
monic potential. Fold. Des. 2(3):173-181, 1997. 

[53] Micheletti, C, Carloni, P., and Maritan, A. Accurate and ef- 
ficient description of protein vibrational dynamics: comparing 
molecular dynamics and Gaussian models. Proteins 55:635- 
645, 2004. 

[54] Zheng, W., Brooks, B. R., and Hummer, G. Protein confor- 
mational transitions explored by mixed elastic network models. 
Proteins 69:43-57, 2007. 

[55] van Kampen, N. G. Stochastic processes in physics and chem- 
istry. 3rd ed., Elsevier, Amsterdam, 2007. 



